Cavity QED in superconducting circuits: susceptibility at elevated temperatures 
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We study the properties of superconducting electrical circuits, realizing cavity QED. In particu- 
lar we explore the limit of strong coupling, low dissipation, and elevated temperatures relevant for 
current and future experiments. We concentrate on the cavity susceptibility as it can be directly 
experimentally addressed, i.e., as the impedance or the reflection coefficient of the cavity. To this 
end we investigate the dissipative Jaynes-Cummings model in the strong coupling regime at high 
temperatures. The dynamics is investigated within the Bloch-Redfield formalism. At low tempera- 
tures, when only the few lowest levels are occupied the susceptibility can be presented as a sum of 
contributions from independent level-to-level transitions. This corresponds to the secular (random 
phase) approximation in the Bloch-Redfield formalism. At temperatures comparable to and higher 
than the oscillator frequency, many transitions become important and a multiple-peak structure 
appears. We show that in this regime the secular approximation breaks down, as soon as the peaks 
start to overlap. In other words, the susceptibility is no longer a sum of contributions from indepen- 
dent transitions. We treat the dynamics of the system numerically by exact diagonalization of the 
Hamiltonian of the qubit plus up to 200 states of the oscillator. We compare the results obtained 
with and without the secular approximation and find a qualitative discrepancy already at moderate 
temperatures. 

PACS numbers: 85.25.-j, 42.50.Pq, 85.25.Cp, 03.67.Lx 



INTRODUCTION 

There is, currently, a substantial activity in the re- 
search into thephysics of Josephson qubits. In particular 
many groups [l|,l2>|3|,|J] became interested in studying the 
systems composed of Josephson qubits and harmonic os- 
cillators, e.g., microwave cavities, mechanical resonators 
etc. Thus the field starts to resemble quantum optics 
where atoms in cavities have been investigated for many 
years. On one hand there are many results which can 
be simply "translated" from the "language" of quantum 
optics to the "language" of solid state physics. On the 
other hand there are specific properties of the solid state 
devices that might require further research. 

In this paper we study the dynamics of a two-level sys- 
tem (qubit) and a cavity at resonance. In quantum optics 
this regime is the most studied and interesting one. The 
corresponding Jaynes-Cummings model (for a review see 
Ref. |5() has been widely investigated in the literature. 
Experimentally, however, the strong coupling regime is 
difficult to achieve and it is also a challenge to keep the 
atom in the cavity for a long time. In optical cavities 
the strong coupling regime was achieved only a decade 
ago 0- Rydberg atoms in superconducting cavities Q 
provide the strong coupling regime and one can even per- 
form quantum gates during the time of flight of the atom 
through the cavity. 

In solid state devices, e.g., a Josephson qubit reso- 
nantly coupled to a damped strip-line superconducting 
cavity 0, the "atom" is permanently placed in the cav- 
ity. This should simplify the time constraints. Also, the 
strong coupling limit seems to be possible. Thus these 
circuits constitute very promising setups for exploring the 
strong coupling limit of cavity QED. Compared to opti- 



cal cavities, one of the differences is the finite tempera- 
ture of solid state devices. Moreover, even if the nominal 
temperature of the refrigerator is much lower than the 
cavity's energy splitting, the hot photons can arrive via 
the leads connecting to the controlling circuits. Thus 
the elevated temperature regime is of substantial inter- 
est. We are aware of only one article addressing the finite 
temperature case, namely Ref. 0- In this paper we ex- 
pand the domain of parameters as compared with Ref. 0. 
Namely we consider the case when both the qubit 's and 
the cavity's dissipative rates are much smaller than the 
cavity-qubit coupling (Rabi-frequency) and both the cav- 
ity and the qubit are coupled to the finite temperature 
baths. Moreover we focus on the correlation functions of 
the cavity, e.g., cavity susceptibility. In solid state sys- 
tems this quantity is particularly convenient to measure 
as it is related to the impedance of the strip-line cavity. 
The temperature-dependence of this impedance consti- 
tutes a direct tool for probing the number of photons in 
the cavity. 

At low temperatures the dynamics is very simple. Only 
a few states are occupied and the susceptibility may be 
presented as a sum of Lorentzians, corresponding to a 
few transitions allowed from these states. E.g., at zero 
temperature only two transitions are relevant and the res- 
onant peak of the uncoupled cavity is Rabi split by the 
qubit to two peaks. At high temperatures the frequencies 
of different allowed transitions are densely packed near 
the oscillator's frequency. This situation is called Liou- 
villian degeneracy as, formally, different modes of the Li- 
ouville evolution operator are almost degenerate |IJ. We 
show that the secular approximation, widely used within 
the Bloch-Redfield formalism, fails due to the Liouvillian 
degeneracy. The insufficiency of the secular approxima- 
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tion was already noticed in, e.g., Ref. In this situation 
one has to take more elements of the Redfield tensor into 
account than required by the secular approximation. On 
the other hand the optical master equation takes all the 
necessary elements into account and thus produces cor- 
rect results. 



EXPERIMENTAL MOTIVATION 

We follow the recent proposal, presented in Ref. 0, 
in which a superconducting charge qubit, coupled ca- 
pacitatively to a cavity formed in a coplanar waveg- 
uide, is shown to be a favorable system for reaching the 
strong coupling regime of cavity QED. Typical param- 
eters are cavity resonance/atom transition frequency of 
luq/2tt — 10 GHz, a vacuum Rabi-frequency of <?/27r = 
100 MHz, a cavity lifetime of 1/k = 160 ns (quality factor 
Q = ujq/k = 10 4 ), and an atom lifetime I/7 = 2/is. The 
system is measured by detecting absorption and phase 
shift of microwaves sent through the waveguide. In this 
paper we calculate the linear response absorption of the 
system. Motivated by the experimental parameters we 
focus on the regime of strong coupling (k,7 <C g), and 
also assume that cavity dissipation dominates over atom 
dissipation (7 < k). We note that the photon temper- 
ature in the cavity may be much higher than the base 
temperature of the cryostat, due to coupling to room 
temperature sources through the waveguide. By measur- 
ing the absorption spectra one may determine the photon 
temperature of the cavity. 



THE JAYNES-CUMMINGS MODEL 

Assuming the qubit to be at the degeneracy point 
0, lll| - and for the moment neglecting dissipation, we 
arrive at the Jaynes-Cummings model, described by the 
following Hamiltonian 

Hjc = -^ q b<J z +awa t a+ | er x (a + a 1 *) , (1) 

where a operates on atom/qubit while a,a^ are ladder 
operators of the cavity/oscillator. We consider the res- 
onant regime when lu osc = uj q b = ujq. The system's 
spectrum for g <C loq is shown in Fig. ^ It is ob- 
tained by, first, analyzing the spectrum for g = (left 
side) and, then, lifting the degeneracies by the coupling 
g (right side). Assuming j « wo we can take into ac- 
count the coupling term only when it lifts the degen- 
eracies in the spectrum. The ground state of the un- 
coupled system is \g) = ||,0), which is non-degenerate. 
All other states are doubly degenerate, i.e. the states 
|t, n) and 1) form degenerate doublets for each 

n > 1. Every such doublet is split by the interaction. 
We can thus define the "bonding" and the "anti-bonding" 
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FIG. 1: Spectrum of the Jaynes-Cummings model at reso- 
nance. 



states \b/a,n) = 2~ 1/2 (||, n) ± \[, n - 1)). The energies 
of these states are E b / an = naj q= y/ng/2. 

OSCILLATOR/CAVITY IMPEDANCE 

In ordinary cavity QED the atom susceptibility is usu- 
ally probed. For the qubit-oscillator system the suscepti- 
bility of the oscillator (cavity) is the most easily measured 
quantity. It is defined as 



X (t) = t9(t)([q(t), q (0)}_) 



(2) 



where q 



v + a. The Fourier transform is x( w ) = 
x'(u) + ix"(ll>) = J dtx(t)e iuJt . The imaginary part of 
the susceptibility is proportional to the dissipative real 
part of the impedance, i.e., x"^) °c KeZ(uj). For a 
closed system with quantum levels \i), with stationary 
occupation probabilities pi, it can be presented as 

x (5(lj — LJfi) — 5((j + Ufi)) , (3) 

where w/,; = (Ef — E^/Ti. Thus it is a series of delta- 
peaks, corresponding to the allowed transitions. 

At T = only the ground state is occupied and the two 
allowed transitions have frequencies uj g - b,i — u$—g/2 and 
^>g\ a,i = ^0 + .9/2 and the matrix elements | (g\ q \b, 1) | = 
I {g\ q \a, 1) I = 1/V2- Thus the uncoupled oscillator's 
peak at ui = ojq is split as shown in Fig. [21 This is called 
Rabi-splitting. 

At T > higher energy states are occupied. The 
transitions are allowed only between neighboring dou- 
blets. There are two classes of allowed transitions. The 
"bonding-bonding" or "anti-bonding-anti-bonding" tran- 
sitions correspond to the transition frequencies 



Ub,n-l; b,n = ^0 ~ ^ (V™ ~ V n - l) , 



^a,n— 1: a,n 



+ 2 (V" _ Vn-l) 



(4) 
(5) 
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FIG. 2: Transitions (left) resulting in a Rabi splitting of the 
cavity susceptibility (right), at T — 0. 
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FIG. 3: Transitions (left) and cavity susceptibility (right) 
at T > 0. The central (blue) arrows correspond to "bonding- 
bonding" and "anti-bonding-anti-bonding" transitions, giving 
rise to the central two peaks in the spectrum. The outer (blue- 
green) arrows correspond to "bonding-anti-bonding" transi- 
tions, giving the outer two peaks in the spectrum. 



with the matrix elements 



KM - 1|<7|M) I = I (a,n- l\q\a,n) | = 



= - (Vn + Vn - 1) 



(6) 



They form the first class of transitions which all 
are positioned inside of the zero temperature Rabi- 
splitting. Analogously, the " bonding-anti-bonding" tran- 
sitions with frequencies 



^>a,n-l; b,n — <^0 

and matrix elements 



W6, n -l;o,n = ^0 + + Vfl— l) , 



+ Vn - 1) 



(7) 

(8) 



| (b, n - 1| q \a, n) \ = | (a, n - 1| q \b, n) \ 
1 



= -(VH-V^I) 



(9) 



are all positioned outside the Rabi-splitting. Both classes 
of transitions are shown in Fig. [3] 



a, n+l^> 

b, n+l> 

a. n > 

b, n > 

a, n-l> 
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FIG. 4: At high temperatures T>uo the "bonding" (blue) 
and "anti-bonding" (red) states form two almost equidistant 
spectrums. The non-harmonicity decreases as \/n + 1 — \fn ss 
1/2-y/n for large n. 



is due to the reduced matrix element for each transition 
as compared to the original oscillator's transition). 

The problem is getting the line-shape at finite (rel- 
atively high) temperatures and high values of Q. In- 
deed, when temperature becomes of order loq more tran- 
sitions become available. To estimate the heights of 
the peaks we note that the occupation probabilities of 
the doublets are estimated as p n ~ Zexp(— uloq/T). 
The matrix elements for the " bonding-bonding" or " anti- 
bonding-anti-bonding" transitions grow as oc \/n, while 
those for "bonding-anti-bonding" transitions decay as 
oc l/y 7 ™. Thus we obtain a "Poisson" distribution for 
the "bonding-bonding" or " anti-bonding- anti-bonding" 
peaks' heights with the maximum at n rs T/ujq. The 
heights of the "bonding-anti-bonding" peaks decay with 
n. Therefore the highest peaks are situated at lu ~ 



ujq ± y'ujQ/Tg/A. 



The spacing between the dominant peaks (~ I 



T I 



thus decrease with temperature, while their width 

increase with temperature. Around the cross-over 
wo(gQ/4wo) 2//5 the peaks start to over- 



temperature T c 
lap. 

Thus, as the temperature grows, we have a transition 
from a quite " un- harmonic" spectrum without Liouvil- 
lian symmetry, where all transition frequencies are dif- 
ferent, to a spectrum which resembles that of two un- 
coupled linear oscillators, see Fig. 0] The widths of the 
peaks in these two cases behave qualitatively different. In 
the first (un-harmonic) case the widths grow with tem- 
perature, while in the second (harmonic oscillator) case 
they do not. 



CAVITY DISSIPATION 

Due to the dissipation the delta functions should be 
widened to Lorentzians. At zero temperature the widths 
of the peaks are estimated as Su> = u>q/(2Q) (the factor 2 



QUBIT /ATOM DISSIPATION 

The qubit/atom is also subjected to dissipation, which 
will add to the peak widths in the susceptibility. At 
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the degeneracy point of the qubit the longitudinal noise 
(coupled to a z ) is suppressed to first order [Tlj. The 
transverse noise (coupled to <r x ) will induce transitions in 
the qubit, characterized by the rate I/7. In terms of the 
Jaynes-Cummings eigenstates the qubit dissipation cause 
similar transitions as the cavity dissipation. A major 
difference is that the matrix elements for these transitions 
are independent of the oscillator state n. Thus, as long 
as the two sources of dissipation are of the same order 
of magnitude, the high temperature susceptibility will be 
determined by the cavity dissipation. 

BLOCH-REDFIELD FORMALISM 



imaginary part of Lx is responsible for the energy shifts 
(Lamb shift). Ly(u>) of the qubit bath is defined and 
treated analogously. 

Considering a high quality cavity (Q 3> 1), and low 
qubit dissipation (7 <C ojq), we limit ourselves to single 
photon exchange with the baths. The possible transi- 
tion energies in Eq. arc u> nn i = ±ujq + 0(g), (see 
Eqs. 0)-©)- Furthermore we assume that the baths 
have no structure on the scale of the qubit-oscillator cou- 
pling g, leaving only one relevant parameter for relax- 
ation \XI Q ) = «/(l - e--/ r -), (Fj ) - 7/(1 - e'^), 
and one for excitation (X 2 _ Uq ) = e -"°/ T " (JQ, (F 2 ^) = 
e-uo/Tyfyz '■ ), for each bath. 



We model the cavity/atom dissipation by coupling an 
observable X/Y of a thermal bath to the g/tr^-coordinatc 
of the cavity/atom: 



H = Hjc + qX + <j x Y + H, 



-Hi 



bath 



bath 



(10) 



One popular choice of a bath is the harmonic oscilla- 
tor one 0, 0. In this case ifbath = an d 
X = J2i x i( b i + b l)- However, as we are assuming weak 
dissipation and doing the lowest order calculation, the 
precise nature of the bath is unimportant and only the 
bath correlator (X(t)X(0)) matters. 

The Bloch-Redfield equation 0, 0] is a kinetic equa- 
tion for the reduced (qubit+cavity) density matrix: 

Pmn(^) ~\~ ^^mnPrnni,^) — ^ ^ Rmnm'n' Pm'n'(t) , (H) 
m' ,n' 

where the Bloch-Rcdficld tensor is given by 



R 



mnm' n 



+ K> 



n' nmm' 



dmm'Kikkn' + finn'Knkkm' 1 (1^) 

k 



and 



Xn'mnn' = (m'\ q \m) (n\ q \n') L x (uJ nn >) 

+ (m'| cr x \m) (n| cr x \n') L Y {u n n>) , (13) 

where Lx (w) is the Laplace transform of the cavity bath 
correlator 



L x (oj) = / dte- lult (X{t)X{Q)) 



(14) 



Introducing the Fourier image of the unsymmetrized cor- 
relator {Xl) = j dt(X(t)X) e iut , we obtain 



L x (oj) 



dv {Xl) 



2tt v + ui — iO 



(15) 



Thus for the real part of Lx which determines the tran- 
sition rates we obtain KeLx(<jj) = (1/2) (X^__J). The 



The secular approximation 

Looking for the slow dynamics of our system, on the 
time-scale given by the dissipation (1/k, I/7), we rewrite 
eq. Hll|l in the rotating frame 



E 



R 



mnm' n 



'(*)■ 



(16) 



where p mn {t) — Pmn{t)e'' UJmrit evolves slowly in time. 

When there are no Liouvillian degeneracies the phase 
(uj mn — Lo m r n :)t in the above expression rotates with a fre- 
quency of the order of luq 3> K, 7, except when m = m! 
and n = nf, or m = m' and n = n'. The secular ap- 
proximation is a random-phase type of approximation, 
keeping only the corresponding elements of Redfield ten- 



sor, R 



ninnin 



and R 



mrann ■ 



Within the secular approximation the Bloch-Redficld 
equation separates into a master equation governing the 
occupation numbers: p nn {t) = R nnmm p m m(t), and a 
simple exponential decay of all off-diagonal elements of 
the density matrix: p mn (t) = p mn (0)e^ l ^ mn+Rmnmn ^ t . 
In the susceptibility we find that the weights of the peaks 
are given by the steady-state occupation numbers, deter- 
mined by the transition rates (Rnnmm)i while the peak 
widths are given by the dephasing rates (Rmnmn)- 



Liouvillian degeneracy 

When difference between two transition frequencies 
u) mn — u> m ',n' ( m 7^ m ' and n ^ n') becomes smaller 
than the transitions' widths oc (k, 7) the peaks start to 
overlap. In this case there is no justification for the sec- 
ular approximation, and more elements of the Redfield 
tensor must be retained. 

In the lowest order (single-photon transitions) the al- 
lowed transition (photon) frequencies are ±wo + 0(g), 
giving two possible values for \uj mn — ^m'n' | = {0, 2u>q} + 
0(g). Looking for slow dynamics we neglect the ele- 
ments Rmnm'n' corresponding to the Liouvillian modes 
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with \u> m n ~ <+>m'n' \ = 2^o + 0(g). This is equivalent to a 
rotating-wave-approximation in the coupling to the bath. 
We keep, however, the elements R m nm'n' corresponding 
to the Liouvillian modes with \w mn — uj m 'n'\ = 0(g). 
That is, grouping the elements of the density matrix 
p mn according to the energy difference E m — E n — 
Mlj{) + 0(g), we now couple only the elements with sim- 
ilar M, but leave the elements with different M uncou- 
pled. We also note that the Bloch-Rcdfield equation with 
this choise of the elements Rmnm'n' is equivalent to the 
"Quantum Optics" master equation: 

p = —i[Hjc,p\ + K.N K (2a)pa~aa)p~paa)) 
+ k(N k + f )(2apa) — a)ap — pa) a) 
+ 7-/V 7 (2(7 + po~ — o~ a + p — p<7~ a + ) 
+ k(N 7 + l)(2cr~per + - (T + a-p - po+ o~) (17) 

where N K / 7 = l/(e u °^ TK ' 1 ~1) are the average occupation 
numbers of each bath at frequency loq. 

NUMERICAL APPROACH 

In our numerical studies we take into account the low- 
est N states of the harmonic oscillator, and then check for 
convergence by increasing N . We work in the eigenbasis 
of the Jaynes-Cummings Hamiltonian (Eq. EJ, which we 
obtain by numerical diagonalization. 

We then rewrite the Bloch-Redfield equation (Eq. Ijlljll 
in a matrix form, i.e. 

P[mn] (t)+iU[ m n][mn]P[mn](t) = R [>nn][m'n'} P[m'n']{t) , 

[m'n'\ 

(18) 

where p(t) — p\ mn ](t) now is a column vector of length 
(2N) 2 and R = R[mn][m'n'] is a matrix of size (2N) 2 x 
(2N) 2 , and uj — ijJ[mn][mn] is a diagonal matrix of the 
same size. The solution can then be written 

pit) = e (--+«)* • p(0), (19) 

which we evaluate by exact diagonalization of — ilu + i?, 
which is the bottleneck of the calculation since the size 
of this matrix grows with A^ 4 . 

Fortunately the property of the Redfield tensor to cou- 
ple only elements of the density matrix p mn with similar 
energy difference E m — E n w Mluq, makes the matrix R 
block-diagonal. The size of each block is only AN x AN , 
which makes the problem tractable up to hundreds of 
states in the cavity. 

Cavity susceptibility 

The cavity susceptibility is defined as 

/>oo 

= i / e*"*<g(t)«(0) - q(Q)q(t))dt, (20) 
Jo 



where the system is in its steady-state at t = 0. In our 
case the steady state density matrix is diagonal with oc- 
cupation numbers, determined by the temperatures of 
the baths, on the diagonal (p s teady)- 

Using the quantum fluctuation regression theorem we 
find 

(q(t)q(0))=Tr{q-fi(t)}, (21) 

where p,(t) is the solution to the Bloch-Redfield equation 
with initial condition /x(0) = q ■ p steady 

Using the solution in Eq. H19J) we can perform the 
Laplace transform in Eq. I|20(l analytically. We now note 
that the main contribution to the susceptibility around 
lu w luq, which is the main topic of this paper, comes 
from the block of R connecting elements p mn with en- 
ergy difference E m — E n sa luq. Thus we only need to 
diagonalize this single block of size 47V x AN to obtain 
the susceptibility. 

NUMERICAL RESULTS 

The numerical results in Figs. l5l()l show the typical be- 
havior of a Rabi-splitted peak at zero temperature, the 
appearance of a broadened multiple peak structure at in- 
termediate temperatures, and the merge of the peaks into 
a single sharp peak at high temperatures. The parame- 
ters used in Fig.|5]are taken from Ref.|3 g = O.OIljo, an d 
Q = 10 4 ,7 = 0.08k. We assume the same temperature 
in both baths. The dash-dotted (red) line correspond to 
the secular approximation. 

If a higher quality cavity can be achieved, more peaks 
will be seen before they merge. This is shown in Fig. 
where Q = 10 5 , and with unchanged qubit dissipation 
(implies 7 = 0.8k). 

CONCLUSIONS 

In this paper we have studied the regime of cavity QED 
relevant for the Josephson qubits in superconducting cav- 
ities. Namely, we considered the strong coupling regime 
when the Rabi splitting is much bigger than the inverse 
life times of the qubit and the cavity, but the temperature 
is high. This regime may be realized due to the hot pho- 
tons penetrating the cavity from the manipulating cir- 
cuits. We have calculated the susceptibility of the cavity 
and found that it is very sensitive to the temperature. 
Thus this quantity might be used for temperature mea- 
surements. On the formal, theoretical side we noticed 
that the secular approximation widely used in the ap- 
plications of the Bloch-Redfield formalism is insufficient 
when Liouvillian degeneracies are present. We compared 
solutions found with and without the secular approxima- 
tion and showed qualitative difference between the two 
at elevated temperatures. 
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FIG. 5: Cavity absorption close to the resonance frequency 
uio, for different temperatures. The paramters are taken from 
Ref. 0, i.e. qubit-cavity coupling g = O.OIljo, cavity quality 
factor Q — 10 4 , and qubit dissipation 7 = 0.08k. 



FIG. 6: As in Fig.[^]but with improved quality factor of the 
cavity Q — 10 (implies 7 = 0.8k). 
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